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ABSTRACT 


In recent years, a large increase in the number of computer 
controlled systems and processes has dictated the use of devices with 
inherently long time delays. The result has been unstable control 
systems. This thesis develops two methods of feedforward compensation 
which cope with this problem. The Pade' method of polynomial 
approximation is used to formulate a number of transfer functions which 
can be used in the two feedforward compensation schemes. To choose 
the best transfer function fora particular problem, a minimum integrel 
Square error analysis technique is developed for use in the frequency 
domain. Computer simulations are used to verify the effectiveness of 


the two compensation schemes. 
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lL, INTEROP VC iio 


In any control system, with lumped linear elements, the output is 
theoretically characterized by a response which starts immediately after 
the application of a transient input. It is well known that this is not 
really true in many cases, since a number of devices may have a 
transportation lag or other delay associated with them. If the delay is 
small, it can usually be ignored. The system is then described by a set 
of linear differential equations and the classical methods of analysis 
can be used to study the behavior of the system, In recent years, the 
large increase in the number of computer controlled systems and 
processes has dictated the use of devices with inherently long time 
delays. Examples of such devices would be a digital computer controlled 
system, where the computation time could be several seconds, or an 
automated superfractionator, which could have a time delay as long as 
several hours. A typical system, with unity feedback, is shown in Pig. I. 

When the effects of time delay cannot be ignored, the primary 
problem is one of system instability. This is easily seen by looking at 
the Bode diagram shown in Fig. 2. This diagram represents the open 


loop transfer function for the system in Fig. 1] with 


| 


S(s+t) 


G, (s) 


and 


G; (Ss) = = ee (2) 
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PLANT WITH TIME DELAY 
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Bode Diagram for Equation 1. 
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The time delay, T, was set equal to 1 second. The upper phase angle 
curve represents the system without the time delay. It can be seen 

that the effect of the delay is to increase the phase lag as the frequency 
increases. The amount of the increase is equal to the product of the 
frequency and the time delay. In this particular example, a very stable 
phase margin of 58.5° has been reduced to 24°. If the time delay is 
increased to 2 seconds, the system would be unstable. It can be easily 
concluded that even a modest delay can present formidable stability 
problems. 

This paper presents a method of compensating for time delay in 
closed loop systems, such as the one shown in Fig. 1. The technique 
involves feeding a conditioned signal forward, around the element 
associated with the time delay. The feedforward path may require either 
active or passive elements to condition the signal. The circuit require- 
ments will be developed. This method is an offspring of one referred to 
yee. 0. BaeLley” ae as the Smith Linear Predictor, However, Buckley's 
use of the predictor was not to stabilize a closed loop system. He was 
concerned with the effects of an extemal disturbance (noise) being 
introduced into the loop between the controller and the plant. He 
assumed the system to be stable before the disturbance was applicd. The 
predictor was only a method for reducing the effect of the distuveame: 
Thé Smith Linear Predictor, as presented in Buckley's paper, is shown 
ig, 3. 

It can be observed that the circuit requirements for the signal 


conditioner in the feedforward path cannot be synthesized with lumped, 
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linear, passive network elements. A suboptimal network will therefore 

have to be used in the feedforward path. The Pade' method of polynomial 
approximation will be derived and used to find an acceptable feedforward 
network. It will be found that a number of acceptable approximations 

can be formed by the Pade' method to meet the network requirements. 

The problem will be to choose the best approximation and to find the 
conditions for which it is the best. In this thesis, a frequency domain 
technique is devised for choosing the best approximation in a minimum inte- 
gral square error sense. Experimental examples are then given to show 


the effectiveness of the approximations and the compensation method. 
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II, FEEDFORWARD COM PENSATION OF 


SYSTEMS WITH TIME DELAY 


Methods of Compensation. Given a system such as the one shown in 
Fig. 1, the solution to the stability problem involves the removal of 

the e Sl term from the denominator of the Closed loop transfer function. 
Iwo methods of doing this immediately come to mind and are Shown in 
Pig. 4. The transfer function for System 1 (Fig. 4a) can be expressed 


as 


COS) GIG SIO" 
Ris) “" 1 + GIG, +GEG west ©) 


It can be seen that if G2 (s) is set equal to Chto! the transfer 


function reduces to 


CO) _ Gs (5) Ga (5) B7Sh 
RES) W i G, (SG Gas} 


It is recognized that this equation represents nothing more than the 
transfer function of a closed loop system (without a delay clement in 
the loop) cascaded with the time delay. 
In system 2 (Fig. 4b), if G(s) iS ag@ingcer Caual to Go(1-e SF), 
the transfer function can be immediately written as ; 
C(s} _ 6, {s) GelS) _ 
Ris} + G(s)Gals) 


It can be seen that both techniques meet the objective of Kelme@ving the 


(S) 





e ST term from the denominator of the closed loop transfer function. 


Both systems use a standard network, which has a transfer function 
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i-e , cascaded with a Gy (s) network. This means that the time 
delay network needs to be synthesized only once, normalized with 
respeci to the time delay, for all problems. The only requirement is an 
impedance match between the time delay network, Go (s) network, and 
any other sections in the feedforward path. There is a subtle difference 
in the two methods of compensation however. This difference involves 
the nature of the output. The two systems can be reduced to the 
equivalent block diagram form Shown in Fig. 5. The output of the first 
system is still delayed by a factor T while the output of the second 
system starts to respond immediately after the application of the input 
signal. First impressions would indicate that the second system is 
better than the first since the time delay has been completely eliminated. 
However, it remains to be seen from experimental data that this is the 
case. It is possible that noise attenuation characteristics, speed of 
response, or other desirable parameter might require. 

use of the first system. A practical consideration might also dictate 
the use of the first system. If, for example, the time delayed plant is 
a chemical process, it would be impossible to "mix" the chemical 
Output with an electrical feedforward signal. On the other hand, the 
electrical feedforward signal could be added to the electrical feedback 
signal coming from some output sensor, as required in the first system. 
The Approximation Requirement, The first step in implementing cither 
of these two systems is to synthesize a network which generates the 


sT 


transfer function l-e °*, which is identical to the transfer function 


required for the Smith Linear Predictor. Since it is a transcendental 


il 





function, some form of approximation must be used if the transfer 
function is to be realized by lumped passive components. Obviously 
distributed parameter delay lines can be used but if Tis large, as can 
be expected in a process control, these devices are noi practical. A 
set of linear filter transfer functions can be derived using any of the 
standard approximation methods. The background, procedure, and 
accuracy behind the Pade' method of polynomial approximation is 
contained in the next chapter. Approximations based on this method are 
used in this thesis for the feedforward network required in the two 
compensation schemes of Fig. 4. 
Summary. The results of this chapter can be summarized as follows: 

1. A closed loop control system with time delay can be stabilized 

if the proper signal can be fed forward around the time delay clement. 

2. The filter required to generate the feedforward signal has the 

transfer function 

G(s) = ites 28 
3. Approximation techniques must be used to synthesize the 


required transfer function with lumped, linear components. 
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Il, DESIGN OF OPTIMAL FEEDFORWARD FILTERS 


The Pade' Approximation. A required step in the synthesis of any 
function, in lumped, linear, passive, time invariant form, is to insure 
that the function is rational. The Pade' approximation is a method of 
representing a polynomial or series of polynomials as a rational 
eo . 
function. Assume that some function, G(s), is to be approximated 
by a rational fraction of an n-th order polynomial in s, divided by an 
m-th order polynomial ins. G(s) could, for example, be the Taylor 


series expansion of the desired transfer function ]-e75!. Thus G(s) 


can be represented, in general rational form, as 


} 
Aj,t QA, St ++ +AnS 
G(s) md ooere: ~ ee aw (6) 


bot b, St °° * + bms™ 





This expression may be expanded by cross multiplication to give 


: (7) 
, a 
= (Agta, StrasrS@¢) 


0 


of G(s) around the origin in the s-plane. By equating equal powers of 


where G iG, etc. are the coefficients of the Taylor series expansion 


s, a set of simultaneous algebraic equations are obtained. Fora given 
number of terms in the Taylor series and for a particular number of terms 
in the rational fraction, there is a finite number of equations which are 


solved simultaneously to obtain the a's and b's. 


is 





From the <0 terms Ag =b.G 


0-0 
1 - Bias : 
2 = 
3 = si 
n = im =e 
S Be ao biGo “I ee her] wy 8 OG + boG, (8) 


The rule for finding the s™ equation should be obvious. It is noted that 
it requires m+n+l terms of the Taylor series to find n+l (orn) coefficients 
for the numerator and m (or m+1) coefficients for the denominator of the 
rational fraction, and that m+n simultaneous equations are used from (8). 
This still allows one of the coefficients, say by, to be picked 

arbitrarily and the only difference in the approximations will be a 
constant multiplier for different values of bo: It is realized that when 
the rational fraction is divided out, in ascending powers of s, the first 
m+n+l terms of the Taylor series are reproduced exactly. Therefore, 

the "maximum errof" incurred in using the Pade’ approximation is no 

more than the remainder of the Taylor series, taken after m+n+l terms. 
Nomenclature, The following nomenclature will be used in the remainder 
of this thesis to refer to a particular approximation or transfer function. 
Peay rerers to the Pade’ approximation, of the function get characterized 
by an n-th order polynomial divided by an m-th order polynomial in s. 


wh 


But e ** is not the transfer function required in the feedforward path of 


the systems in Fig. 4. Therefore, G,,, will refer to the approximation 


~siT 


of the desired transfer function, l-e ~~, derived from the Pade' 
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e ne “of cr * 
approximation of e Sl’ of the same order. For example, if 


Aa 

= 
xy } 
Ow | 
las 


then 
657 +5°T* 
Gy == = = eee (10) 
El e Lice f ee ape & 
6 7 i | spelt } 7, | 
Derivation of the Desired Filter. The Taylor series expansion of e *, 
about the value x = 0, is the well known power series 
—X 2 S 3 (11) 
x = eco Mv ~———— rfooume GOS 
This is the series that is used to derive rational approximations for the 


=S i = 


function l-e , with st =x. Pade’ approximations, for the function e 
are tabulated in many sources.” Appendix l(a) gives such a list of 
approximations. Appendix 1(b) gives the derived filter transfer 
functions derived from the appropriate approximation, while Appendix 1(c) 
gives the derived filter transfer functions in factored form, which is not 
available in the general literature. The polynomials were factored using 
the POLYRT subroutine on the CDC 1604 digital computer and rounded 

off to three decimal places. While the factored form is never used 
directly in any calculations in this paper, it does give a little insight 

as to the location of the poles and zeros of each of the filter transfer 
functions. The poles and zeros are plotted in Fig. 6. There appears to 
be a pattern to the poles and zeros as the order of s in either the 


numerator or denominator is increased. This pattem is not known but 


if it were, it could be used to formulate a general rule for finding the 
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filter transfer function of any desired order. However, the solution 
of (8) is quite straight forward in solving for any particular Ban so that 
no effort was expended to determine an expression for this pole zero 
pattem. 
There is no guarantee that all of the transfer functions found in 
Appendix 1(b) are realizable with passive eclerents. The properties of 
a transfer function, made up of Jumped, linear, passive elements are 
well known but are reiterated here. © 
1. Poles are in the left half of the s plane or, for the lossless 
case, on the imaginary axis and simple. Complex poles must be 
conjugate pairs, 
2. zeros may have any s plane location but complex zeros must 
occur in conjugate pairs. A minimum phase transfer function has 
all of its zeros in the left half of the s plane. 
3. Poles may not be located at the origin or at infinity. 
4, The degree of the numerator may equal but not exceed the 
degree of the denominator. 
Item 4 above rules out all transfer functions where m is less than n. 
These transfer functions will no longer be considered in this thesis. 
Referring back to Fig. 6, it is seen that not every Pade' derived transfer 
function can be synthesized with lumped, linear, passive elements, 
even though m is greater thann. As the order of s in the denominator 
is increased, holding the order of s in the numerator constant, the 
motion of the poles is such that poles begin to appear in the right half 
of the s plane. This violates rule 1 above and renders these transfer 
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functions unacceptable as an approximation to the desired filter 
transter funccioen, 1-e75! It is also noted that for every transfer 
function where m equals n, an amplifier is required, for D.C. gain, in 
addition to the passive elements for synthesis. 

There is no indication as to which of the filter transfer functions 
Po. ithe best.” Intuitively, Gjz7 should be the best eppre moi ienmrosac 
ideal transfer function since its derivation required more terms of the 
power series than any of the others in the list. On the other hand, Ga 
requires more effort to synthesize since it is more complex than the | 
others. It may be that a simple filter will work just as well as a 
complex one in the feedforward path of the systems in Fig. 4. A 
measurement of which is the best approximation to the ideal filter 


sT 


transfer function, l-e °*, is attempted below. 


The Minimum Integral Square Error Index. Tiewind ex, ecnesen stommeasuxe 


moaen of the filter transfer functions is the best approximation, is 


the Minimum Integral Square Value of the difference between the 


sT 


responses of the ideal transfer function, l-e °°, and the Pade' derived 


filter approximation G_,,. The integral square error = be calculated in 
either the time or frequency domains. Although the time domain has 
received predominant attention in the literature, several advantages 
uncovered during this investigation urge the use of the frequency domain 
technique. The first, and most important of the advantages, is the 

ability to solve for a frequency spectrum which can be later used to find the 
integral square error for any input. This spectrum, called the weighting 


function, gives the user a certain feeling for the "quality" of the 
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network under consideration. When time domain techniques are used, 
an input signal must immediately be chosen and the result applies for 
that input only. The result that is obtained is only a number and is 
not a very satisfying "answer." Without other networks (numbers) for 
comparison, the results are next to useless. Even with several networks, 
the user only knows which one has the smallest integral square error, not 
what causes that error or what effect another input might have. This 
advantage will be more apparent by comparison of the results of the 
two techniques. Another advantage of the frequency domain approach 
(or disadvantage of the time domain method) is not having to take the 
inverse LaPlace transformation of the network transfer function. Some 
of the larger Pade' approximations are very difficult and time consuming 
when it comes to obtaining time domain solutions. 

In view of the above arguments, this thesis will primarily use the 
frequency domain approach to selecting the best Pade' approximation 
for use as a filter transfer function in the compensation systems of 
Fig. 4. It is known that the mean square value of a signal is equal to 


the value of its autocorrelation function at t = 0, y 


mean sq. value of v(t) = Om (o) 
CP ' 
=] (TUT He) dT 
La 
0 oo 
= ia Oyy (7) SU 
= a 


where & (iw) is referred to as the Power Density Spectrum of the 
signal. It derives this name from the fact that the integral in (12) equals 
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the average power dissipated in a 1 ohm resistor by the signal 
voltage v(t). @ yy iw) can be interpreted as a lineardensity distribution 
of the mean square value of the signal along the imaginary axis of the 


ST 


s plane. (),,,(jw) is obtained from the relationship 


>< oS, 
Ow a Ase - fe Oy (C} IT 02) 


em (%D 
and then sctting s = jw. 


Consider the hypothetical experiment shown in Fig. 7. It can be 
shown’ that the power density spectrum of the error signal Y(s) can be 


expressed as 


Oyy(s) = H(S)HES) O, (5 
- HEIG(S)6.2(9 
- His) Gf-s) O,,(s) 
+ G6(s)G(-S) Qs) (14) 


where 
.. @ § 
Hfs) = le € a 


and O11) is 1/Ziitimes the LaPlace transform of the autocorrelation 
function of v,(t ), 6 12(s) is 1/Zitimes the LaPlace transform of the cross 
correlation function of vj (t) ‘and Poll, Gwen 1 eyl(e) EUS Sap liyl, Se 
shown by the dotted lines in Fig. 7, then all of the auto and cross 
correlation functions in (14) are equal and will be designated by a 
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Therefore, from (14) with s = jw 


Oy, Guts Oy Gal lute “+ ley * 
“ez Ra dH (2upE" (pu) }}- y 


The quantity in brackets in (17) will be referred to as W(jw) 


WisW)= Hiiwl “e 6g Se Relive ida 


By trigonometric substitution 


off 


H(gm=zl-@ 0 sl fees (a wt) td sively 7) 


s(l-cos WT) + ¢ SiN wT 


so that 


. 12 2 
| H (gu) = @(|—cos wi 


From the Pade’ approximations, G 


as) is of the form 
A= 
Gl ) _ fn S4 2 Ay) 9 | + ES ENG # Ao 


— wo f 
Bin SB, Braet S s° cect BS 4- Ee 


where m and n depend on the particular approximation being 


considered. Orin simpler terms 


Ms) 4 - Ny, (S s) 
Me (s) } Ne (S} 


ZS 


G(s) 





(16) 


(17a) 


(18) 


(19) 


(20) 


(21) 





where the M's refer to the sum of the terms with even powers of s and 
the N's refer to the sum of the terms with odd powers of s. For real 
frequencies (when s = jw) 


| Mle+F Mh MY | 
Gye) = pone = Dtwl+é Ble? som 
MoT pis (us) 


and 


Fr 


joreuy A= + Be (23) 





From (18) and (22), it follows that 


Re[ Maw) 6*(gey)] 
= A(w)l- COs WT] BAL sew wT] 


Substituting (19), (23) and (24) into the bracketed portion of (17) yields 


(24) 


Wise) =2 (i- cos wT) hie ty (25) 
2 [Aly (i-cos wt) +B0Nsiv Wi] 


In shortened form, (17) becomes 


OD yy (Zu) = Wfew) Oy (eu) (26) 


W(jw) is given the name Spectral Weighting Pilea. Referring to 

Fig. 7, the spectral weighting function does just what its name implies. 
It weights the magnitude of the power density spectrum of some input 
signal v(t) to yield the power density spectrum of the error. All that is 
now required, to determine which of the Pade' derived filter transfer 
functions is the best, is the substitution of the various approximations, 


Gmn- into (25) to determine W(jw). Multiplying the spectral weighting 
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function by the power density spectrum of an appropriate input signal 
yields the power density spectrum of the error for that particular 
transfer and input function. Integrating the power density spectrum 
according tO (12) yields the integral ee value of the error function. 
The Pade’ derived filter transfer function that has the smallest area 
under the spectrum, i.e. the minimum integral square error, is the one 
which is the best approximation to the ideal transfer function, Pere 
The Spectral Weighting Function. In order to generalize the selection 
process for any input signal, the spectral weighting functions, for each 
mime acceptable Pade derived transter functions , Should be presenmed 
in graphical form. This permits the user to know in advance which 
frequencies of the input signal the various transfer functions weight 

the most. The transfer function, that best suppresses the frequencies 

in the input signal, is the one that will have the minimum integral square 
error. It is significant to note that the best approximation depends 

upon the spectral distribution of the input signal as well as the spectral 
weighting function. Due to the large number of Pade’ derived filter 
transfer functions and the complexity of the spectral weighting function 
equation, these curves were computed using a digital computer. A 
program was written to compute and plot W{jw) and is presented in 
Appendix 2(b). The results are shown in Appendix 3. The horizontal 
axis for the graphs in Appendix 3 is plotted as the logy 9 of a normalized 
frequency time delay product, i.e. x = 10919 wT. In general, the 
spectral weighting functions for the Pade' derived filter transfer 
functions have a larger weigiting effect at higher frequencies and 


ot 





appear to act like high pass filters of cutoff frequency, w With 


0: 
the exception of the curves corresponding to the case when m equals n, 
all of the spectral weighting functions settle to a constant amplitude of 
one. The spectral weighting function of a transfer function having m 

equal ton has a constant amplitude oscillation about the value two. 

If the frequency components of an input signal were lower than some Wg) 
where the spectral weighing function is zero, then the integral square 
error would be zero. For such an input signal, the Pade' eae filter 
transfer functions would be an excellent approximation to the ideal 
transfer function. 

Comparison of Filter Transfer Functions. The input signal that is quite 
often used to test the performance of a control system or filter is the 

unit step. The idea is that any signal can be approximated by a sum 

of little steps. Any filter or system that would perform well with a 

single step input should work for a sum of steps. A unit step may not, 
however, be a very good choice for an input signal in the selection 
process described above. Note that the minimum integral square error is 
a nonlinear performance index and the principle of superposition cannot 

be used. The Pade' digenes filter transfer function that is selected as 

the best, based on a unit step, cannot be assumed to be the best for 

some other input. The minimum integral square error performance index 
does not, however, rule out the possibility that the best transfer 

function for a step input may also be best for other inputs. It all depends 
on the power density spectrum of the input. For the purpose of this paper, 
a unit step will be used as the input signal when selecting a Pade' 
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approximation to use as the feedforward filter transfer function. This 
choice was made simply because the unit step is general and casy to 
work with. In an actual problem, where the input signal is known, ithe 
power density spectrum of the actual input would be used. Also, the 
unit step makes the examination of experimental results a familiar task. 


The power density spectrum of a unit step is given by the equation 


0, Mm 4 Wy] = wager (27) 


‘the digital computer is also used to multiply the spectral weighting 
function by the power density spectrum of a unit step and to numerically 
integrate the resulting area under the power density spectrum of the 
error, The program that was written to perform this operation is given in 
appendix 2(c). The integral square error for each of the acceptable Pade' 
derived Pisin di functions is tabulated in table 1. The table is broken 
into several bands to show how the frequencies of the spectrum cont- 
ribute to the integral square error, Figs. 8 and 9 are graphical 
representations of the data in table 1. The solid lines in Fig. 8 show the 
variation in the integral square error as one of the indices (m or n) is 
varied. Horizontal movement along one of the solid lines from one 

point to another represents a change in the order of s in the numerator, 
while holding the order of s in the denominator constant. It is observed 
that, as long as nis less than m, an increase in the order of s in the 
numerator results in a decrease in the integral square error. The amount 
of the decrease becomes smaller as n approaches m. When n becomes 


equal to m, there is a substantial increase in the integral square error and 
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when n becomes larger than m the integral square crror goes to infinity. 
Moving vertically from one solid line to another represents a change 

in the order of s in the denominator while holding the order of s in the 
femerator constant. For these changes, the #meegral squarewemremaveus lia 
decreases as the order of s is increased, The exceptions to this 

occur when there is a pole "close" to the imaginary axis (Fig. 6) such 

as G4o and Gel . The dotted lines in Fig. 8 represent the change in the 
integral square error as both indicies are changed, but the difference 

is walle constant. It is readily observed that the amount of improvement 

in the integral square error is exponentially decreasing as the order (and 
complexity) of the transfer function is increasing. Figs. 9 (a and b) 

show the changes in the contribution to the integral square error by the 
various frequency bands. It is observed, from Fig. 9(a), that changing 
the order of = in the denominator, while holding the order of s in the 


numerator constant, has little effect on which band contributes the most 


M.to the integral square error. On the other hand, Fig. 9(b) shows 


that changing n while holding m constant has a considerable effect on 
the frequency band contributing the most area. As nis increased, the 
higher frequency bands tend to contribute more. This fact could be used 
to advantage in many systems. A large number of closed loop systems 
tend to act as a low pass filter. If one of the Pade’ derived filters is 
used in such a closed loop system, the low pass filtering effect of the 
system should suppress the higher frequencies that are weighted more 

by the high pass feedforward filter. This means that the integral square 
error for a unit step input to such a system would be less than that 
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when n becomes larger than m the integral square error goes to infinity. 
Moving vertically from one solid line to another represents a change 

in the order of S$ in the denominator while holding the order of s in the 
numerator constant. For these changes, the integral square error usually 
decreases as the order of s is increased, The exceptions to this 

occur when there is a pole "close" to the imaginary axis (Fig. 6) such 

as G4o and G6 1° The dotted lines in Fig. 8 represent the change in the 
integral square error as both indicies are changed, but the difference 

is held constant. It is readily observed that the amount of improvement 
fmeene Integral square error is exponentially decreasing as the order (and 
complexity) of the transfer function is increasing. Figs. 9 (a and b) 

show the changes in the contribution to the integral square error by the 
various frequency bands. It is observed, from Fig. 9(a), that changing 
the order of s in the denominator, while holding the order of s in the 


numerator constant, has little effect on which band contributes the most 


Mereo the integral : square error, On the other hand, Fig. 9(b) shows 


that changing n while holding m constant has a considerable effect on 
the frequency band contributing the most area. As nis increased, the 
higher frequency bands tend to contribute more. This fact could be used 
to advantage in many systems. A large number of closed loop systems 
tend to act as a low pass filter. If one of the Pade' derived filters is 
used in such a closed loop system, the low pass filtering effect of the 
system should suppress the higher frequencies that are weighted more 

by the high pass feedforward filter. This means that the integral square 
error for a unit step input to such a system would be less than that 
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calculated above. It would therefore be advantageous to use the filter 


transfer function with the highest frequency weighting in systems that 


act as low pass filters. 


fhe Integral Square Error Index ,(Time Domain). In order to check the 


above results, the integral square error is calculated for a few cases 


in the time domain. The same hypothetical experiment, shown in Fig. /7, 


is used. The equation for the error signal is 


cm & 5P 
4). = 
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Taking the inverse LaPlace Transform 
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The integral square error is given by the equation 
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It is easily seen that the first integral reduces to T (the area under a 
pulse of magnitude 1 and duration T). The second integral is equal to 
the area under the response curve of the Pade' derived transfer function, 


G from time Oto T. The last integral represents the area under the 


mn’ 
square of the filter response curve for all time. The lower limits can 
be changed to 0 since the response does not start prior to the application 


of the input at time 0. The equation for the integral square error, for a 


unit step input, therefore reduces to 


(33) 





The digital computer program presented in appendix 2(d) was 
written to compute I from (33). Because of the long computation time, 
only a few of the Pade' derived filter transfer functions were tested by 
this method, The results are given in table 2. Since all of the 
approximations contain decaying exponential terms, the upper limit for 
the third integral in (33) was set at 20 times the longest time constant. 
By comparing tables 1 and 2, it is seen that the results of the two 
methods agree to within the third decimal place. Better agreement can 
be obtained if the tolerance on the numerical integration is made smaller, 
Such action would increase the computation time however. 

Other Approximations for y-e 7S! Several other sets of approximations 
were investigated but none performed as well as the Pade'. The defining 


equations for the other approximations investigated are® 
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7 (j- As)! 
aad fp S) i (35) 
These were chosen mainly on “bce that they are relatively easy to 
synthesize. For example, (35) is a lattice network yielding the proper 
time delay. The results of the investigation are given in table 3. The 
format is the same as that for table 1 (i.e. the frequencies in each 
band are the same as forthe bands in table 1). It is noted that the integ- 
ral square error for the Bo approximation becomes smaller as nis 
increased. However, for the same degree of complexity, the Pade’ 
approximation has a smaller mean square error, The C, approximation 
does not perform well at all. The smallest integral square error occurs 
when n= I and as nis increased, the integral square error increases. The 
smallest integral square error TO approximation is larger than the 
largest integral square error for a Pade! approximation. 
Summary. From the above observations, it appears that the Pade’ 
approximation derived from the largest number of terms in the power 
series has the smallest integral square error and is therefore the best 


sl . This statement 


approximation to the ideal transfer function, l~-e~ 
is, of course, subject.to the condition that the denominator is of higher 
order than the numerator and that there are no poles in the right half of 
the s plane. It is speculated however, that the use of this approximation 
may not be the most practical in many cases, The increase in complexity 


(and therefore in difficulty in synthesis) may far outweigh the advantage 


of a small decrease in the integral square error. As stated before, the 
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selection of how complex a Pade' derived filter transfer function to use 
will depend on the system being compensated, the duration of the time 
delay, the nature of the input signal, etc. The results of this chapter 


can be summarized as follows: 


1. The Pade’ approximation is a method for representing a 
transcendental function as a rational fraction of polynomials. It 
is found by the simultaneous solution of a series of algebraic 
equations. 

2. Some of the Pade' approximations of the desired function 1] eros 
can be realized with lumped, linear passive elements. 

3. The best approach in choosing a filter transfer function is* 

through the frequency domain. In this procedure, a quantity called 
the Spectral Weighting Function is generated from the desired transfer 
function and the transfer function of a filter under consideration, 

The spectral weighting function graphically shows which frequencies 
in the input signal contribute most to the integral square error of the 
filter. When the spectral weighting function is multiplicd by the 
power density spectrum of an input signal, the power density 
spectrum of the error is obtained. The filtcr with the minimum area 
under the power density spectrum cf the error is the one with the 
minimum integral square error and is considered "the best." 


4, The Pade' approximations have the smallest integral square error 


of all of the approximations investigated for the step input. 
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IV. EXPERIMENTAL INVESTIGATIONS 


In order to examine the effects of approximation in the feedforward 
compensation schemes outlined earlier, a hybrid computer model was 
tested with several of the Pade' derived filters in the feedforward path. 
The specific objectives of the tests were to ascertain the derogatory 
effects of using an approximation, rather than the ideal transfer function 
ee, and to determine what effects the order of the approximation 
has on the time domain performance of the feedforward filter. 

The Plant. For simplicity, the plant described in the introduction was 
the one used for the hybrid computer model, The plant is shown again 
in Fig. 10(a). Fig. 10(b) shows the computer model. An EAC model 
TR-20 analog computer was used to simulate the linear portions of the 
plant. A CDC model 160 digital computer was used as the plant time 
delay element. In addition, a portion of the analog computer provided 
signal conditioning for proper A/D and D/A conversion. The A/D 
converter requires all signal inputs to be between 0 and -10 volts. The 
analog computer setup is shown in Fig. 11. The maximum allowable 
Signal at any point in the computer model is chosen to be +8 volts. 
This is a convenient choice of limits since the digital computer 
Operates on the octal number system and a direct correlation between 
voltage levels in the analog computer and numbers in the digital 
computer memory results. 


The Digital Time Delay. The basic idea behind the digital time delay 





is not complex. The analog computer Output is sampled periodically 
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and the number corresponding to the amplitude of that sample is stored 
in the computer memory. At the end of the time delay Same T, the 
number is returned to the analog computer as a voltage signal through 
a hold element. If the sampling rate is fast enough, the woitsee 
returned to the analog computer appears to be a continuous time varying 
signal, A flow graph of this operation is shown in Fig. 12. The 
computer program written to perform this operation is given in appendix 2(f). 
There is a difficulty in the use of the program listed in appendix 2. (f) 
however. As the time delay is increased, the sampling rate is decreased, 
This program is only adequate for for delays of several minutes or less, 
For time delays longer than this, the output of the D/A converter is better 
described as piecewise constant rather than continuous. 
Simulation of a Typical Pade' Derived Filter. The most difficult step in 
setting up a computer model of the system is the simulation of the 
feedforward filter, Cee Large gains (larger than 107 in Go 4) are 
required in all but the simplest of the Pade' tee filters. The steps 
taken to simulate Gy3 are given here as an example of the procedure. 
__ 81057 + 6O(STIRLAO (ST}*e (sth 
G43 ~ pun augg ST OSTe. Pals 7) °4-1G(ST 7)Bts 7} 


Dividing the numerator and denominator of (36) by (st) 4 and setting 


(36) 


equal to the ratio of the output to the input. 
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By cross multiplication and rearranging terms 





(38) is in the proper form for the simulation on the analog computer. 
The analog computer setup for G42 is shown in Fig. 13. The other Pade' 
derived filter approximations are simulated in the same way. 
The Step Response. It was stated in the introduction that the system 
under consideration was stable with a time delay of 1 second but 
unstable if the time delay was increased to 2 seconds. The Bode 
diagram for the plant was presented in Fig. 2 to support this. To verify 
this, the feedforward filter was disconnected in the computer model and 
the model was subjected toa 1 volt step. The results are shown in 
Fig. 14. The plant is definitely unstable for a time delay of 2 seconds, 
To meet the objectives of the tests, three of the Pade' derived 
filter transfer functions were used as the feedforward network in the 


computer model. The three that were used are COE G43, and Gece 


These were chosen because they had the minimum integral square error 
for a given value of m. The time delay was set to 2, 3, and 4 seconds. 
At all three of the time delay settings, the plant without compensation 
would be unstable. The results of the tests, for both cornpensation 
methods, are Sheu in Figs. 15 (a, bandc). It is observed that all 
three Pade’ derived filter transfer functions stabilize the system for all 


three values of time delay tested. For the longer delays, the more 
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complex filters do a better job with respect to the shape of the response 


curve. The response with the G,-. filter and with a 4 second time delay 


Oo 
looks just like the response of the uncompensated system with no time 
delay, with the exception of the time shift. It is interesting to note that 
the two compensation methods (Fig. 4) do not give the same compensated 
system output. The best example of observing this is the response of 
Goy with a 4 second delay, Fig. 15({a). The response due to compensa- 
tion method 2 oscillates much more than the response due to compensa~ 
tion method 1 and the settling time is longer. On the whole, it appears 
that compensation method 1 performs better than compensation method 2, 
This is very fortunate since it was noted earlier that practical 
considerations might require the use of this compensation scheme. 
Sensitivity. No extensive investigation was made into the area of 
sensitivity. Only one short experiment was performed to get an idea of 
the tolerance limits on the filter time delay setting. The plant was 
simulated on the hybrid computer with a two second time delay. This 

is an unstable condition. The plant was then compensated, using 
compensation method 1] and the Pade' derived filter Goy° The filter time 
delay was varied from 1 to 5 seconds, in half second steps. The 1 
second setting was the shortest delay that could be conveniently set 

on the computer and at 5 seconds, the system was at the margin of 
stability. The results of some of the other trials are shown in Fig. 16. 
It is observed that the system is stable for a large range of settings of 
filter time olay The quality of the response deteriorates as the 
difference between the plant and filter time delays increases but the 
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V. CONCLUSION 


It was easily shown, in Chapter 2 of this thesis, that the 
feedforward compensation techniques presented in Fig. 4would remove 
the phase lag due to a time delay in a closed loop control system and 
would therefore improve the stability of such a system. It was then 
implied that the required feedforward network could not be synthesized 
with lumped, linear, passive, time invarient circuit elements. The 
problem was therefore reduced to one of finding a rational approximation 
that would perform as near to the unrealizable ideal circuit as possible. 
Since the function to be synthesized is a common one, a number of 
approximations and approximation techniques were available to choose 
from. The performance index, chosen to rate the various approximations, 
was the minimum integral square error between the step response of the 
ideal and approximation filters. The Pade' method of polynomial 
approximation _ finally used since it was simple to apply and it gave 
a wide variation in transfer function complexity. This large variety of 
circuits means that a large range of integral square error was available to 
test the effects of approximation in the two compensation schemes, It 
was found, through computer simulations, that approximations of the 
ideal transfer function effectively stabilize systems with time delay. 
Even poor approximations, i.e. filters with large integral square errors, 
were able to effect a degree of stability in such systems. In addition, 
it was found that the stabilizing properties of the two compensation 


schemes were not effected by small changes in the time delay setting 
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system remains stable. This feature is useful in situations where the 

plant time delay is either not known exactly or it changes. 

Summary. The observations made in this chapter can be summarized 

as follows: 
1. The Pade' derived filter transfer functions make effective 
feedforward filters in the compensation schemes outlined in Fig. 4. 
As anticipated, the more complex transfer functions, which are 
better approximations to the ideal transfer function, give a better 
system response. 
2. Compensation method 1 (Fig. 4(a)) gives a better system response 
than compensation method 2. 
3. The ability of the Pade’ derived filter transfer functions to 
stabilize an unstable plant is not very sensitive to changes in the 


filter time delay setting. 


oi 





V, CONCLUSION 


It was easily shown, in Chapter 2 of this thesis, that the 
feedforward compensation techniques presented in Fig. 4would remove 
the phase lag due to a time delay in a closed loop control system and 
would therefore improve the stability of such a system. It was then 
implied that the required feedforward network could not be synthesized 
with lumped, linear, passive, time invarient circuit elements. The 
problem was therefore reduced to one of finding a rational approximation 
that would perform as near to the unrealizable ideal circuit as possible. 
Since the function to be synthesized is a common one, a number of 
approximations and approximation techniques were available to choose 
from. The performance index, chosen to rate the various approximations, 
was the minimum integral square error between the step response of the 
ideal and approximation filters. The Pade' method of polynomial 
approximation was finally used since it was simple to apply and it gave 
a wide variation in transfer function complexity. This large variety of 
circuits means that a large range of integral square error was available to 
test the effects of approximation in the two compensation schemes. It 
was found, through computer simulations, that approximations of the 
ideal transfer function effectively stabilize systems with time delay. 
Even poor approximations, i.e. filters with large integral square errors, 
were able to effect a degrce of stability in such systems. In addition, 
it was found that the stabilizing properties of the two compensation 


schemes were not effected by small changes in the time delay setting 


is) 





of the approximation feedforward filter. In general, it can be concluded 
that both the compensation methods presented in Chapter 2 are excellent 
compensators for time delay and that the required feedforward circuits 
are not very sensitive to approximation effects. 

The concept of a frequency domain analysis of the minimum integral 
square error criterion, as developed in Chapter 3, is as important as 
the idea of feedforward compensation. In most of the literature, the 
time domain analysis is used. While the time domain technique gives 
satisfactory answers, i.e. gives the integral square error for a problem 
under consideration, the frequency domain analysis gives this 
information and much more, The speciral weighting function, an inter- 
mediate step in the frequency domain analysis, gives the user a 
graphical indication of the frequencies that are contributing to integral 
square error, This means that the user can find a solution for his 
problem based on the frequency spectrum of the input signal and the 
spectral Pe sciina functions of the various solutions under consideration, 
The spectral weighting function can be computed once and used to 
calculate the integral square error for any input. These advantages are 
not available when the time domain analysis is used. 
Topics for Further Study, The ideas used in this thesis can be projected 
into additional areas of study. For example, work can be done in the 
case where the time delay cannot be isolated from the remainder of the 
open loop transfer function of the plant. Also, the idea of a sensitivity 
investigation could receive additional attention. It is interesting to note 


that the ideal feedforward transfer function is very similar to the 


3), 





transfer function of a Zero Order Hold circuit, (l-e7S!)/s. A number 

of the Pade' derived filters found in appendix 1 could possibly be used 
in conjunction with an integrator to produce a passive approximation of 
a zero order hold circuit. The nadaiencierm ues are a promising method 
for finding a rational approximation for any irrational function and work 
can be done in the implementation of such functions. Another area 
where work can be done is to consider phase relationships in the 
frequency domain analysis of the error. Only the amplitude relationship 
(integral square error) was considered in this thesis. It would be 
interesting to see what additional information could be obtained from 


such an analysis. 
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